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Abstract 

By studying the effective potential of the MSSM at finite temperature, we find that CP 
can be spontaneously broken in the intermediate region between the symmetric and broken 
phases separated by the bubble wall created at the phase transition. This type of CP 
violation is necessary to have a bubble wall profile connecting CP conserving vacua, while 
violating CP halfway and generating sufficient baryon number without contradiction to 
the experimantal bounds on CP violations. Several conditions on the parameters in the 
MSSM are found for CP to be broken in this manner. 
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1 Introduction 



The idea of the electroweak baryogenesis [p]] is attractive in that it could solve the problem 
of matter-antimatter asymmetry in the universe by the knowledge of accessible exper- 
iments on the earth. In particular, the nature of CP violation, which is one of the 
requirements to generate the baryon asymmetry of the universe (BAU) starting from the 
baryon-symmetric universe, will be revealed in the near-future experiments. The viable 
mechanisms of the electroweak baryogenesis, however, require some extension of the stan- 
dard model with other sources of CP violation than the phase in the Kobayashi-Maskawa 
matrix. 

Among the extensions, the minimal supersymmetric standard model (MSSM) may not 
only admit various CP violations but also cause first-order electroweak phase transition 
(EWPT) with the small soft-supersymmetry-breaking parameter in the stop mass-squared 
matrix]^, |3|. It is also pointed out that the chargino and stop may play important roles in 
transporting the hypercharge into the symmetric phase, where it biases the sphaleron pro- 
cess to generate the BAU[|]]. The CP- violating phases in the mass matrices are essential 
in these scenarios, while they are constrained by various observations such as the neutron 
electric dipole moment (EDM). Another source of CP violation, which was originally con- 
sidered in the baryogenesis mechanism ||, is that in the Higgs sector, that is, the relative 
phase of the expectation values of the two Higgs doublets. The Higgs VEVs including the 
phases, which characterize the expanding bubble wall created at the first-order EWPT, 
vary spatially. This spatially varying CP violation makes, through the Yukawa coupling, 
the quarks and leptons to carry the hypercharge into the symmetric phase 0. This sce- 
nario will work even if the superpartners are so heavy that they are not excited in the 
hot plasma. Of course, since this CP-violating phase also enters the mass matrices of the 
charginos, neutralinos, squarks and sleptons, this might enhance the generated baryon 
number when they are thermally excited to act as the charge carriers. 

In a previous paper |7|] we attempted to determine the profile of the bubble wall by 
solving the equations of motion for the effective potential at the transition temperature 
(> Tc — lOOGeV) in the two-Higgs-doublet model. For some set of parameters, we 
presented a solution such that CP- violating phase spontaneously generated becomes as 
large as 0(1) around the wall while it completely vanishes in the broken and symmetric 
phase limits. We shall refer to this mechanism as 'transitional CP violation'. This solution 
gives a significant hypercharge flux, by the quark or lepton transport ||. We also showed 
that a tiny explicit CP violation, which is consistent with the present bound on the 
neutron EDM, does nonperturbatively resolve the degeneracy between the CP-conjugate 
pair of the bubbles to leave a sufficient BAU after the EWPT|J. 
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In this paper, we examine the possibility of the transitional CP violation at finite 
temperature in the MSSM by calculating the effective parameters, which are defined as 
the derivatives of the effective potential. Similar analyses were executed by use of the 
high-temperature expansion of the finite-temperature corrections, one of which concerned 



the problem of the order of the EWPT[12] and another focused on the spontaneous CP 



violation in the broken phase[|TTJ. It should be noted that the high-temperature expansion 
is not always a good approximation, especially when the masses of the particles running 
through the loops are larger than the transition temperature. We apply it only to the 
light stop loop, while the contributions from the other particles are treated numerically. 
Although a tiny explicit CP violation is needed to have nonzero BAU||, we shall con- 
centrate on the possibility of spontaneous CP violation. In § 2, we briefly review the 
mechanism of the transitional CP violation. We derive the formulas for the effective pa- 
rameters, which include both zero- and finite-temperature corrections, in § 3. We show 
the numerical results and analyze the possibility for the spontaneous CP violation in § 4. 
Discussions are given in § 5. The calculations of the loop corrections and the relevant 
integral formulas are summarized in Appendix. 



2 Transitional CP violation 



(2.1) 



Consider a model with two Higgs doublets whose VEVs are parameterized as 

= (! = 1 ' 2) 

and 9 = 9\ — 9 2 . We assume that the gauge-invariant effective potential near the transition 
temperature has the form ofQ 

1 



' 2 2 , 1 2 2 2 /) i ^1 4 , "^2 4 

2 m iPi + 2 m 2^2 - m 3 p lP2 cos 9 + — Pl + — p 2 



+ 



A3 + A4 2 2 A 5 2 2 1 2 2 

"Pl/°2 + -rPlP2 cos2 " ~ ni^Pl + A 7 p 2 )pip 2 COSfc' 



4 ^r. 4 

[Ap\ + PiP 2 (A) + B x cos 9 + B 2 cos 29) 

+Pipi(C + Ci cos 9 + C 2 cos 29) + Dpi 

TTPiP2 - 2(5 2 PiP2 + C2P1P2) 



X 



cos# — 



2ml + A 6P? + A 7 p| + 2(B lPl + C x p 2 ) 



2A 5 pip 2 - 8(B 2P i + C 2 p 2 ) 
-l-^-independent terms. 



(2.2) 



1 All the parameters in V c g should be regarded as the effective ones containing both zero- and finite- 
temperature corrections. 
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The p 3 -terms are expected to be induced at finite temperature in a model whose EWPT 
is of first order. Since we do not consider any explicit CP violation, all the parameters 
are assumed to be real. For a given (p 1; p 2 ), the spontaneous CP violation occurs if 

F(pi,to) = ^ P 2 lP l-2(B 2 plp 2 + C 2Pl pl)>0, (2.3) 
,^ rfn n x 2m l + Vi + A 7 p| + 2(B lPl + C lP2 ) 

-l<Lr{ P l, P2 = — — — — — r < 1. [2 A) 

2A 5 pip 2 - 8{B 2 pi + C 2p2 ) 

In the MSSM at the tree level, A 5 = A 6 = A 7 = and A = B k = C k = D = (k = 0, 1, 2), 
so that no spontaneous CP violation occurs. At zero temperature (A = B k = C k = D = 
0), it is argued that A 5 6j7 are induced radiatively and (|2.3| ) is satisfied if the contributions 
from the chargino and neutralino are large enough. For (|2.4j) to be satisfied, m\ should 
be as small as \% P \ + A 7 p| so that the pseudoscalar becomes too light fll3| , 



At T ~ Tc, the values of (pi, p 2 ) vary from (0, 0) to (vc sin /3c, vc cos (3c) between the 
symmetric and broken phase regions, where the subscript C denotes the quantities at the 
transition temperature. Then the effective parameters in (|2.2|) include the temperature 
corrections as well. Hence there arises large possibility to satisfy both (|2.3|) and (|2.4j) in 
the intermediate region at the transition temperature, without accompanying too light 
scalar. If this is the local minimum or a valley of V e g appears for intermediate 

{puPz) with a nontrivial 9. It should be noted that such a local minimum need not 
to be the global minimum of the effective potential. For such a V e g with appropriate 
effective parameters, the equations of motion for the Higgs fields predict that some class 
of solutions exist, which have 9 of 0(1) in the intermediate region even if it vanishes in 
the broken phase p. 

In the following sections, we calculate the effective parameters in (|2.2j ) to examine 



whether the conditions ( j2.3|) and (2.4) are satisfied for some intermediate ( P i, P 2)- 



3 Effective parameters of the MSSM 

Since we are concerned with the possibility of the spontaneous CP violation, all the 
parameters in the lagrangian are assumed to be real. The tree-level Higgs potential of the 
MSSM is 



V = m x ip\ip d + m 2 (fl(p u + (m 3 (f u (p d + h.c) 



+ 



(3.1) 
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where 



ml 



Ai 
A 5 



m 2 d + \p\ 2 , 
^2 = \(g 2 2 + 9 2 i), 



2 ~ 2 i I |2 

m l = m u + \p\ 



mn 



A3 = -Xg\-g\ 



pB, 
A 4 = 



Ar — A? 



0, 



(3.2) 
(3.3) 



Here g 2 (i) is the S77(2)([/(l)) gauge coupling, /i is the coefficient of the Higgs quadratic 
interaction in the superpotential. The mass squared parameters m 2 ud and pB come from 
the soft-supersymmetry-breaking terms so that they are arbitrary at this level. m\ could 
be complex but its phase can be eliminated by the redefinition of the fields when A 5 = 
A6 = A7 = 0. We adopt the convention in which this m\ is real and positive. 
Let us parameterize the VEVs of the Higgs doublets as 



J_[Pi 
V2\0 



— ( 



-( ° 



The effective potential at the one-loop level is 



v eS = v + v 1 (p i ,e) + v 1 ( Pi ,e;T), 

where Vi(pi,0) is the zero-temperature correction given by 



(3.4) 



(3.5) 



m. 



64tt 2 



log 



ML 



(3.6) 



and Vi(pi,9;T) is the finite temperature correction; 

T 4 r°° r / / m 

V^pi, 0;T) = —Y / rij J q dxx 2 log [l - sgn(n,) exp (-^x 2 + m 2 /T 2 )J . (3.7) 

Here we used the DR-scheme to renormalize V e s with the renormalization scale M rcn . 
rij counts the degrees of freedom of each species including its statistics, that is, rij > 
(rij < 0) for bosons (fermions). mj, which is a function of the Higgs background (pi,0), 
is the mass eigenvalue of each species. 

At the one-loop level, {m\) eE receives corrections only from the Higgs bosons, squarks, 
sleptons, and charginos and neutralinos. A 5) 6,7, which are zero at the tree level, are gener- 
ated only through the loops of these particles. Among them, we consider the contributions 
of charginos(x ± ), neutralinos(x°), stops(t) and Higgs(0 ± ). The effective parameters are 
defined as the derivatives of V c s at the origin of the order-parameter space: 

d 2 V of 



(ml) 



/ eflf 



eff 



dvidv 2 



m 3 + A x m 3 + A^m 3 + A^±m 3 , 



(3.8) 
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A 5 = -r 
A 6 = 

A 7 = 



<9 4 K 



eflf 



2 \ dv\dvl 
1 9 4 Kff 



d 4 V r 



eff 



dv\dvl 



A X A 5 + A t -A 5 + A^±A 5 , 



3 dv\dv2 
1 <9 4 \/ cff 



3 <%i<9t> 2 



= A X A 6 + A f A 6 + A^±A 6 , 



A X A 7 + A f A 7 + A^±A 7 . 



(3.9) 
(3.10) 
(3.11) 



The explicit forms of the corrections in term of the Feynman integrals are summarized 
in Appendix. In the following, we present the formulas for these corrections from each 
species. 

(i) chargino and neutralino 

The mass matrices of the charginos and neutralinos are given by 



M x o = 



M 2 

M 2 


-\92pi 



-p j 


Mi 




/i 



2 9ip2e 
H 




-id 



(3.12) 
(3.13) 



/ 



respectively. As noted in Appendix, all the contributions from the neutralinos are propor- 
tional to those from the charginos, when the gaugino mass parameters satisfy M 2 — Mi. 
We assume this for simplicity. Then the corrections from the charginos and neutralinos 
are summarized as follows: 

1 



A x m t 



A X X 5 



9 2 2 



\ COS z Vw / 



1 + 



8tt 2 



1 + 



A X A 6 = 



9$ 



8tt 2 



1 + 



cos 4 9 W 
2 



7T 2 T 4 



1 + 



cos 2 9w 
2 



,( +) /M 2 /i\ 

(3.14) 



cos 4 9 W 



+ 



9_l 

7T 2 



cos 4 6 W 
2 



(3.15) 



1 + 



COS 4 



pM 2 {+) /M 2 //\ // 3 M 2 (+) /M 2 /i 



V T ' T 



y4 



V T ' T 



(3.16) 



where L(M 2 ,/i), if (a), ii(a;) and / 2; 3^(a,&) are defined in Appendix, 
(ii) charged Higgs bosons 

The mass-squared matrix of the charged Higgs bosons is 



Mh 



9 2 2 -9l n 2 , 92+91 n 2 



m{ + ^Pt + 
. m l + \9lp1P2e 



-P2 



m l + \9lpip2e ie 
ml + ^pj + ^pl 



(3.17) 
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The low-energy parameters in this matrix are arranged to break the gauge symmetry, 
so that one of the mass-squared eigenvalues evaluated at pi = should be negative; 
m\m\ — m\ < 0. This negative mass squared makes the finite-temperature corrections 
to the effective potential complex for small p^. (Suppose negative m 2 in (|3.7|) .) This 



pathology will be cured by taking the higher-order corrections into account |14[ . Among the 
corrections, the so-called 'daisy diagrams' are the most dominant ones at high temperature 
since they grow as T 2 . Hence we replace m\ and m\ in the Higgs loops with the 'daisy- 
corrected' ones given by 



m^ 



m\ + k{T 2 



m n 



mi + k 2 T z 



(3.18) 



where 



fci 

k-2 



1 



167T 2 
1 



(3s 2 2 + gl 



16n 



2 (3g 2 2 +gl + 4y*) 



(3.19) 



This determines the limiting temperature T\ ow , under which the origin of the effective 
potential is not a local minimum, that is, m\m\ < {m 2 ) 2 : 



rp2 

low 



{k\m\ — k 2 m 2 ) 2 + Akik 2 {m 2 ) 2 — {k\m\ + k 2 ml) 
2k t k 2 



(3.20) 



This limiting temperature is rather large for tan/?o > 2. It will be shown numerically 
that the Higgs contributions to the effective parameters are smaller by order four or five 
than the others as long as they are well defined. Hence we simply neglect the Higgs 
contributions when T < T\ ow in the following calculations. Even if the approximation 
by simply substituting m\ 2 for m\ 2 in the Higgs loops is not justified, the origin of Kfr 
should be a local minimum at T ~ Tc as long as the EWPT is of first order. Then the 
effective m\m\ would be large enough so that the contributions from the Higgs would 
become small, as is clear from the following formulas. 



A (f) ±m 3 
A0±A 5 

A,/,± Ar 



-g 2 m 3 L(p u p 2 ) + -^g 2 m 3 f. 



2_2/(-)// i l ^2 



T T 



92 m t R ( ftf 
64vr 2 plpl \pl, 



4 4,(-)// 1 l ^2 



8vr 2 T 4 



4 2 
92 m 3 




-H 



A, 



+ 



2p? cos 2 9 



T T 
1 



(3.21) 
(3.22) 



8tt 2 T 2 



1 



2 cos 2 8 W 



^2 



K 



A, 



A-)(V± /fa \ fjh _ 

/ 3 y T , T J^{ T2 2T2cos 2 6w 



1 - 



2 cos 2 6 W J T 2 



mt 



)ft ] 



T ' T 



(3.23) 
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8vr 2 T 2 



J 6 



T ' T 



Vy2 



2 cos 2 flu// T 2 



m 2 . 



2T 2 cos 2 6 W 



(-) 



(Hi 


H2\ 




T) 



where 



Ml,2 



m 2 + ± J (ml — ml) 2 + 4m| 



(3.24) 



(3.25) 



The definitions of the various functions used in the above formulas are given in Appendix, 
(iii) light stop and p 3 -term 

When m 2 in fl3.7|) vanishes as pi — > 0, the second and higher derivatives of it for the 
bosonic loops are ill-defined at pi = 0. This singularity originates from the zero mode in 
the summation over the Matsubara modes. Upon approximated by the high-temperature 
expansion [Q, (|3.7|) receives p 3 -terms with positive coefficients from the bosonic particles 
whose masses behave as 0(p 2 ) for pi ~ 0. This p 3 -terms are supposed to make the 
EWPT first order. In the MSSM, the candidates generating such terms are the weak 
gauge bosons, the Higgs bosons and the scalar partner of the quarks and leptons with 
appropriate mass parameters. Among them the Higgs bosons and the squarks and sleptons 
could yield ^-dependent p 3 -terms, that is, B\ i and/or Ci t a in (2^2), which will affect the 
conditions ( |2.3[ ) and ( |2.4| ), if their mass eigenvalues vanishes as pi — > 0. When one of 



the soft-supersymmetry-breaking mass parameters in the squark mass matrix vanishes, 
this situation is realized. Here we consider only the top squark (stop) because of its large 
Yukawa coupling. 

The mass-squared matrix of the stop is 



M 



m 
m 



where 



m 



ii 




92 j O1-P2) + \vtPl, 



m 



-22 



rn 



2 1 2/2 2 \ 1 22 

m t + q9i(Pi- P 2 ) + 2 ytp21 
Vt 



\:2 



V2 



[(ppi + A t p 2 cos 9) - iA t p 2 sin 1 



(3.26) 



(3.27) 

(3.28) 
(3.29) 



Here m|, m 2 and A t come from the soft-supersymmetry-breaking terms and y t is the top 
Yukawa coupling. Although the relative phase between /1 and A t yields an explicit CP 
violation, we assume they are real. 
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The temperature-dependent part of the stop contribution to the effective potential is 



V~ t {p u 6; T) = 3^ \2I B {a 2 + ) + 2I B (a 2 



2tt 



(3.30) 



where Ib{.& 2 ) is defined by 



I B (a 2 ) 



J dx x 2 log (l — e" 



(3.31) 



and the factor 2 counts the degrees of freedom of complex scalars and a± = m±/T 2 with 



m i 



m ll + m 22 ± V i m ll ~ m 22) 2 + Mm 



121 



(3.32) 



being the eigenvalues of M 2 . Since /i and A t are real, \iA t ^ is required to give 6- 
dependence in m 2 ± . If mq and/or m- t vanishes, m 2 , and/or m 2 _ behave as 0(p 2 ) for p ~ 0. 
We assume m 4 ~ ~ and 3> T, since too light left-handed stop, which couples to the 
SU(2) gauge bosons, would lead to a large correction to the p-parameter. Then m± is 
given approximately as 



m 



m 



m \ + o ( -V + d I (pi ~ Pi) + T^Vtpl + 



XY Y 2 



Y 



^m 2 - Arriq 8m. ~ t 



2 i 2 
1 2/ 2 2 \ i 2 2 

gPi(Pi -P 2 ) + 2^^2 - 



XY Y z 



Y 



4m? 4m~ 8m|^ 



(3.33) 



(3.34) 



where 



<?i+2 2 (P1-P2), 



y = 2y 2 (p 2 pj + 2pA t p 1 p 2 cos6 + Afp 2 ). 



(3.35) 
(3.36) 



To evaluate the stop contributions to the effective parameters defined in ( |3.8| ) - ( |3.11| ), 
we need the behavior of V[(pi,8;T) at pi ~ 0. For this purpose, we employ the Taylor 
expansion around a~ = mhjT 2 to evaluate Ib(o? + )\ 



IbK) 



^B( a l) / 2 \2 



+ ..., 



(3.37) 



where a 2 + = m 2 + /T 2 . On the other hand, we use the high-temperature expansion for 



4 2 

1 a„ a„ + A_a„ + . . . 

45 12 p " 6 p " p ~ 



(3.38) 
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where up to 0(a 3 ) the coefficients are obtained from the well-known formula [[14]], which 
includes a 4 log a-term in addition to a 4 -term. In exchange for dropping the a 4 log a-term, 
we have decided A_ = 0.1764974 by numerical fitting for < a_ < 1. 
We obtain the stop corrections to the effective parameters: 



Arm 2 



A t ~A 5 



N c yffiA t L( m<l ,0)+ , 



3T 2 y^A 



16tt 2 m\M m 



ArA 



i A 6 



ArA 



t A 7 




(3.39) 



(3.40) 



(3.41) 



A- ,(3.42) 



where N c = 3 and Mm is the infrared cutoff parameter, which will be taken to be the order 
of the transition temperature This is needed because of the infrared singularity encoun- 
tered in the presence of a massless particle through the loops, as is well known |i~5|j . This 
is cured by calculating the fourth derivatives away from the origin. Then, by minimizing 
the effective potential, the mass scale in the logarithm is replaced by the VEV, that is, 
the dimensional transmutation occurs. We have checked that as long as Mm > lOOGeV, 
M/^-dependence is not so significant that we simply use Mm instead of minimizing the 
V eS . 

Now we are ready to extract p 3 -term in the stop contribution, which is given by. 



Vi(Pi,0;T) 



0(p3) 



rp4 

3^ I B (a 2 



O(aS) 



rp4 



3/2 



(3.43) 
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When At/rriq, n/rriq, {gi/yt) 2 *C 1, we pick up higher-order terms to obtain 



T 



0(p3) 



+ 



\Vt\ 



1 - 



9{ 



3Mt 



3/i 2 ^ 
2mf, 



P1P2 + 

i(A) 

2 I m J 



3/V 

2 \mqy 

3 (fiA t 



+ 



6 \yt. 



p\p 2 cos 2 9 
2" 



p x p\ cos0 + ---U3.44) 



From this expansion, we extract 



Ci 



B 5 



T 



4^7T 

T 
4V2tt 



■ 3pA 
[3 / //>l t 



1 f A 



6 \vt. 



2 \ m„- 2 



(3.45) 
(3.46) 



Because 



pA t jmh <C 1, it is expected that B 2 is much smaller than C\. We can neglect 
B\ and C 2 , which will be induced from the sbottom loops, compared with C\ and B 2 , 
because the bottom Yukawa coupling \y b \ is much smaller than \y t \. 



4 Numerical Results 

We examine whether the conditions G(pi, p 2 ) < 1 an d |G(pi,p 2 )| < 1 are satisfied or 
not by evaluating the effective parameters included in F and G. As note above, we 
numerically calculates the integrals in the finite-temperature corrections. Before showing 
the numerical results, we comment on some general properties of the behavior of the 
parameters. 

If only the light stop contributes to the ^-dependent p 3 -terms, C 2 = so that F(pi, p 2 ) > 
holds for p 2 satisfying 

(A 5 p 2 - AB 2 ) Pl > 0. (4.1) 
As long as we take p\ to be positive, this implies 

ABo 

P2 > -tA for A 5 > (4.2) 
^5 

ABo 

— - < p 2 < 0, for A 5 < 0. (4.3) 
A5 

The latter case corresponds to a negative tan j3 at finite temperature. Since we adopt 
tan (3q > at the tree level in the following examples and we found several CP-violating 
bubble wall solutions for A5 > 0f|, we concentrate on the former case here. Note that 
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positive corrections to A5 come from the charginos and neutralinos at zero temperature 
and the light stop at finite temperature. All the other contributions are always negative. 



bution is much smaller than those from the charginos and neutralinos. The maximum is 
realized around p 2 = Mf, which corresponds to the maximum of K(M 2 /p 2 ). Since K(a) 
slowly varies around the peak, A x As(T = 0) is positive for a rather wide range of M 2 /p 2 . 
We have checked that even if the finite-temperature corrections are taken into account, 
A 5 is positive for 1/2 < M 2 /p 2 < 3 for T < lOOGeV. In fact, A 5 > 10~ 4 at T ~ lOOGeV 
and B 2 = O(10~ 5 )T ~ lCT 3 GeV in our examples, so that the condition F(p\,p2) > is 
satisfied for p 2 > lOGeV. This will not impose a strong constraint as long as tan /3c is not 
so small. On the other hand, for \G(pi, p 2 )\ < 1 to be satisfied, the value of the tree-level 
777J must be tuned since the magnitude of {jn\) cE must be the same order as \p\ + A7P 2 
and C1P2, which are radiatively induced at finite temperature. 

Now we examine the condition \G(pi, p 2 ) | < 1 for several sets of the tree-level param- 
eters, which are m 2 , m\, m|, A t , M 2 , p, rrig and m 4 ~ = 0, at temperatures of O(100GeV). 
Instead of giving m\ and m 2 , we input the values of the tree- level tan /3q and the absolute 
value of the Higgs VEV t>o, which are related to the masses-squared parameters by the 
relations defining the minimum of Vq\ 



As seen from the formulas for the effective parameters, the signs of the corrections depend 
on those of pM 2 and pA t . For example, the sign of A x m 2 (T = 0) is the same as that of 
pM 2 , while the temperature corrections to it have the opposite sign. If we take \A t \ <C 
\M 2 \ ~ rrig, the chargino and neutralino contributions dominate over those from the stops. 
As long as we adopt a positive m§ at the tree-level, negative pM 2 is needed to have a nearly 
zero {m 2 ) efi at T ~ Tc- Since f 2 T \mi/T,m 2 /T) is a positive and increasing function of 
T, the finite-temperature part of A x m 2 (T) works to reduce {rn\) cS to almost zero from 
its positive zero-temperature value. 

Hence we take two parameter sets with pA t > and pA t < 0, respectively. For 
each case, the T-dependences of the effective parameters are studied and \G(pi,p 2 )\ is 
plotted in (pi, p 2 )-p\a,ne at several temperatures with tan/5 = 1.2 and tan/3 = 5. All the 
numerical values having mass dimension should be understood to be in the unit of GeV. 
We use vq = 246, m t = 177 and M ren = Mjr = 100 in these examples. 



For the expansion ( |3.44|) to be valid, pA t /m^ <C 1 so that we expect that the stop contri- 




m 2 tan /3o - -m 2 z cos(2/5 ), 
m 2 cot (3q + -m 2 z cos(2/9 )- 




(4.4) 



(i) pA t > 

The parameters in the first example is given in Table 0. 
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Ill 

Table 1: The parameters used in the numerical analysis in the case of pA t > 0. 
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Figure 1: {jn\) eS and A 5 as functions of temperature T. The total values are given by 
the solid curves, the corrections from the stop, chargino-neutralino and the charged Higgs 
bosons are depicted by the dashed, dotted and dotted-dashed curves, respectively. For 
T < T low = 71.4 the Higgs contributions are ignored. 

(1) tan A, = 1.2 

For these parameters, T[ ow = 71.411, below which we neglect the contributions of the 
charged Higgs bosons. As shown in Figs. [I] and [|, their values are much smaller compared 
to the others so that we expect them not to alter the results significantly. In this case, 
we have 

B 2 /T = 1.532346 x 10" 5 , C x /T = -4.846485 x 10" 3 . (4.5) 

As seen from the curves in Fig. ||, A 6 7 = O(10~ 3 ) so that C\p 2 is comparable to A 6 pf + A7P2 
for pi ~ 100 at T ~ 100. Hence, when |G(pi,p2)| < 1, ( m i) e fr can ^ e l ar g er compared 
to the case of spontaneous CP violation at T — if C\ has the same sign as it, which 
is the present case. From ( |3.45|) , C\ has the opposite sign to pA t . At T = 100, A5 = 
6.11026351 x 10~ 4 so that F{p u p 2 ) > for p 2 > 4B 2 /X 5 = 10.03. \G{p u p 2 )\ is plotted 
for T = 70, 75, 80 and 85. There exists a region where the condition \G(pi,p 2 )\ < 1 is 
satisfied, as shown in Figs. |3| and [| at each temperature. 
(2) tan (3 = 5 

Now Ti ow = 312.48, so that we ignore the charged Higgs contributions. The behaviors 
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Figure 2: A6 and A7 as functions of temperature. 



T=70 IG(p r p 2 )l 



200 




T=75 

200- 



IG(p r p 2 )l 



150 



100 



50- 



0- 




1.8 



1.6 



1.4 



1.2 



50 100 150 200 50 100 150 200 

Figure 3: Contour plots of |G(pi,p 2 )| at T = 70 and 75. |G(pi,p 2 )| is satisfied in 
black region. 
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Figure 4: The same as Fig. [3] at T = 80 and 85. 
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Figure 5: {inl) cS and A5 as functions of temperature. The total values are given by the 
solid curves, the corrections from the stop and chargino-neutralino are depicted by the 
dashed and dotted curves, respectively. 



0.0005- 
0- 

-0.0005- 
-0.001 - 

-0.0015 
-0.002 

-0.0025 
-0.003 



0- 

-0.0005- 
-0.001 - 

-0.0015- 
-0.002 

-0.0025 
-0.003 

-0.0035 



50 



100 



150 



100 



150 



Figure 6: A 6 and A 7 as functions of temperature. 



of the effective parameters are qualitatively the same as the example above, as depicted 
in Figs. [5] and || Tthe contributions from the charginos and neutralinos are identical to 
those above as obvious from ( |3.14j ) - ( |3.16|) . The larger tan/3 implies the smaller y t for 
a fixed m t , that is, the smaller |A^m|(0)|, which implies the larger (ml) cS for pA t > 0. 
This lowers the temperature at which \G(pi, p 2 )\ < 1 is satisfied for some (pi,p2). In this 
case, we have 

B 2 /T = 7.368276 x KT 6 , d/T = -2.313638 x 10~ 3 , (4.6) 

which implies F(pi,p 2 ) > for p 2 > 5.104 at T = 100, that is, F(p 1 ,p 2 ) > holds in the 
whole region in which \G(pi, p 2 )\ < 1 as shown in Figs. [7| and |8|. 
(ii) pA t < 
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Figure 7: Contour plots of \G{pi, p2)\ at T = 80 and 90. \G(pi,p2)\ is satisfied in the 
black region. 
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Figure 8: The same as Fig. [7| at T = 100 and 110. 
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Table 2: The parameters used in the numerical analysis in the case of pA t < 0. 
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Figure 9: (mj) eff and A5 as functions of temperature T. The total values are given by 
the solid curves, the corrections from the stop, chargino-neutralino and the charged Higgs 
bosons are depicted by the dashed, dotted and dotted-dashed curves, respectively. For 
T < T[ ow = 75.3 the Higgs contributions are ignored. 



Although the stop contributions change their signs, those from the charginos and neu- 
tralinos are still dominant for the parameters in Table |2|. This makes the temperature 
dependences of all the effective parameters milder than those in the case of \xA t > 0. That 
is, for a wider range of temperature, the conditions for the spontaneous CP violation will 
be satisfied. 
(1) tan/3 = 1.2 

For this, T\ ow = 75.313, below which the Higgs contributions are neglected. The effective 



parameters are plotted in Figs. and 10. Since 



B 2 /T = 3.447779 x 10" 



d/T = 7.269727 x 10" 



(4.7) 



F{PiiP2) > for p 2 > 16.05 at T = 100. The almost whole region where \G(pi,p 2 )\ < 1 



satisfies this condition as well, as shown in Figs. [11] and [12]. The positive C\ requires 
smaller m| than that for C\ < to make 2 (m|) efl + 2C\p 2 nearly equal to \p\ + \-jp\- 
(2) tan A) = 5 

Now T low = 313.08. We completely ignore the Higgs contributions in the plots of the 
effective parameters in Figs. [0| and [14]. By the same reasoning as in the case of pA t > 0, 
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Figure 10: Xq and A7 as functions of temperature. 
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Figure 11: Contour plots of \G(pi,p2)\ at T = 90 and 100. \G(pi, p2)\ is satisfied in the 
black region. 
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Figure 12: The same as Fig. 11 at T = 110 and 120. 
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Figure 13: (Tnf) eS and A 5 as functions of temperature. The total values are given by the 
solid curves, the corrections from the stop and chargino-neutralino are depicted by the 
dashed and dotted curves, respectively. 




Figure 14: A6 and A7 as functions of temperature. 
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Figure 15: Contour plots of \G(pi, p 2 )\ at T = 60 and 70. p 2 ) | is satisfied in the 

black region. 
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Figure 16: The same as Fig. [15] at T = 80 and 90. 



the temperature at which \G(pi,p2)\ < 1 is lowered. At the same time, A5 becomes larger 
at lower temperatures so that the region with \G(pi,p2)\ < 1 grows as seen from Figs. [15 



and [16[ For this parameter set, we have 

B 2 /T = 1.657862 x 10" 5 , d/T = 3.470457 x 10" 3 , 
which implies F(p u p 2 ) > for p 2 > 8.463 at T = 100. 



(4.t 



5 Discussions 

We have investigated the possibility of the new type of spontaneous CP violation, which 
occurs at finite temperature in the transient region from the symmetric phase to the 
broken phase separated by the electroweak bubble wall. Since this type of CP violation 
disappears in the broken phase at zero temperature, it is free from any constraint on CP 
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violation from the experiments. Further it will enhance the generated baryon number 
by the electroweak baryogenesis mechanism. Although the CP-conjugated pair of the 
bubbles degenerate in their energies as well as their nucleation rates, a tiny explicit CP 
violation consistent with the observation such as the neutron EDM is sufficient to resolve 
the degeneracy and to leave the present BAU|§. 

For this mechanism to work, some constraints are imposed on the parameters in the 
MSSM. First of all, pM 2 < is required to make {in 2 ) eS decrease by the chargino and 
neutralino contributions as noted in the previous section. \p\ ~ \M 2 \ yields a positive 
correction to A5 by these particles up to T ~ lOOGeV. For the expansion used for the 
light stop contributions to be valid, p 2 < m? and A 2 < m| are required. Together 
with \p\ ~ \M 2 \, these conditions make the stop contributions smaller than those from 
the charginos and neutralinos. The ^-dependent p 3 -terms are induced only if one of 
the soft-supersymmetry-breaking masses of the stop almost vanishes; mi ~ 0. Among the 
coefficients of these terms, C\ gives contributions to the numerator of G(pi, p 2 ) comparable 
to those from the charginos and neutralinos. Whether |G(pi, P2) I < 1 is satisfied is 
sensitive to {mf) efi in its numerator. Since 1C\p 2 appears in the numerator with the 
same sign as (m 2 ) off , the effect of positive C\ (pA t < 0) is compensated by reducing the 
tree-level m\ as shown in the two examples in the previous section. 

One might wonder how the difficulty encountered in the case of the spontaneous CP 



violation at zero temperature |13| is avoided. A5 > is satisfied at T = if the chargino 
and neutralino contributions dominate over those from the Higgs bosons and stops. This 
also applies to F(p 1 ,p 2 ) in our case, except for the Higgs boson contribution, which is 
small since the effective m\m\ is larger than m| at T ~ Tc- The problem was that for 
\G(p\, p 2 )\ < 1 to be satisfied at T = 0, where B\$ = C\ )2 = 0, {m\) eS must be so small as 
(A 6 cos 2 P + A 7 sin 2 P)v 2 , which implies that the pseudoscalar boson is too light. Its mass 
rriA is related to m\ at the tree level by m\ = (tan/3 + cot f3o)m\. The smallest value of 
this tree-level tra in our examples is vtia = 67GeV for m| = 2200GeV 2 and tan/3 = 1.2. 
However, tua should be calculated at the minimum of the corrected potential. The value of 
rriA calculated in this way might be sufficiently large for suitable range of the parameters. 
Further the parameters of the theory are not so constrained in our case compared to the 
case at T = 0. This is because, if only the conditions are satisfied at some (pi, p 2 ) in the 
transient region, the CP phase could be large enough to generate the BAU around the 
bubble wall. 

Finally we emphasize that the spontaneously-CP-breaking minimum does not have 
to be the global minimum of V e g. The transitional CP violation could take place if the 
conditions are satisfied for some fixed (pi,p 2 ), since such a bubble wall with transitional 
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CP violation would have a lower energy than that without CP violation. Hence we do 
not need to be afraid that such a local minimum may not be the absolute minimum. 
For another reason, however, we should understand the global structure of the effective 
potential, which determines Tq, to know whether the transitional CP violation occurs or 
not. In this sense, the conditions we examined here should be regarded as the necessary 
conditions but not the sufficient ones. With the knowledge of the global structure of V e e, 
one could find the CP-violating profile of the bubble wall so that one could estimate the 
generated baryon number. 
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A Loop Corrections 

Here we summarize the expressions for the contributions to the effective parameters 
in terms of the finite-temperature Feynman integrals from the charginos, neutralinos, 
charged Higgs bosons and stops at finite temperature. We also show various formulas to 
calculate the Feynman integrals. According to the definitions of the effective parameters 



( p.8|) - (|3.11|) , the contribution of each particle is expressed in terms of the propagators 
in the symmetric phase (pi = 0) and the vertices which are related to the derivatives of 
the mass matrices. 

The contributions from the charginos whose mass matrix is given by ( |3.12| ) to the 
effective parameters are 

(+) 



A x ±m 2 3 = 2glfiM 2 i f A l (k)A 2 (k), (A.l) 

Jk 

A X± A 5 = -2<7 2 4 (//M 2 ) 2 z [ {+) Al(k)Al(k), (A.2) 

Jk 

/•(+) 

A X ±A 6 = A X ±A 7 = 2g\pM 2 i \ k 2 A\(k)A\(k), (A.3) 

J k 



where 



and the integral implies 

_ y?, f d 3 k yfi+Yi J = u*)n = Znin/T for bosons ,^ , 

~ ' (27r) 3 1 ko — icu n = (2n + l)m /T for fermions. 



Ik 



n=— oo 
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The neutralino contributions are rather lengthy because of its mass matrix given by 



A x om 2 



A x o A5 



A x oAg 



(+) 



2(02 + / A 1 (k)A 2 (k)A 3 (k) (k 2 - P 2 )Q 2 + PR 1 

Jk L 

-4Q? 2 2 + g{) V< Al(k)Al(k)Al(k)(k 2 -P 2 - R 2 ) 

Jk 

(k 2 - P 2 )Q 2 + PQR 2] 



x 



A x oA 7 = 2(gj + gl)ni [ W Aj(k) Aj(k) A 2 3 (k)k 2 \k 2 - P 2 - R 2 



X 



{2k 2 - 2P 2 - R 2 )Q + PR 2 



(A.6) 



(A.7) 



(A.8) 



where 



A X (A;) 



^1,2 



2 • 



A 2 (k) 



k 2 — 



2 ■ 



A 3 (k) 



k 2 - jj? 



k 2 -fij 

P 2 + Q 2 + 2R 2 ±\P + Q\ J(P - Q) 2 + AR 2 



(A.9) 
(A.10) 



with 



P = M 2 sin 2 9 W + Mi cos 2 9 W , 
Q = M 2 cos 2 9 W + Mi sin 2 6 W , 
R = (M 2 - M^smOwcosOw. 



(A.ll) 
(A.12) 
(A.13) 



If Mi = M2, P = Q = M 2 and R = so that fii = fi 2 = M 2 . In this case the neutralino 
contributions are reduced to 



A x om 2 



A x oA5 
A x oA 6 



— — — A x ±ml, 
cos 2 6 w x 6 

9 

-A V ±A 



cos 4 fl 



X ± A 5, 



A x oA 7 



- A^±Afi 



-A y ±A 7 . 



cos 4 6 w x cos 4 9 w 

We consider this special case of M 2 = Mi for simplicity. 

The corrections to the effective parameters from the charged Higgs bosons are 



(A. 14) 
(A.15) 
(A.16) 



A^im 2 
A^±A 5 
A^,±A 6 



■^ 2 m 2 z^ ( A 1 (k)A 2 (k), 



4/™2\2 , 



t#2 ( m : 



(-) 

k 

(-) 



A 2 (A;)A 2 (A;), 



(A.17) 
(A.18) 



■-glmjz jf A 2 (A;)A 2 (fc) 
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A^±A 7 



-4 2 • 
j-02™3* 



fc — ( 1 



2 cos 2 #w 
(-) 



2 cos 2 9w 



1 

2 cos 2 9w 



ml 



2 cos 2 9 W 



where 



2 

^1,2 



1 



k 2 /i 2 
mf + m\ ± 



m i) 2 + 



(A.20) 

(A.21) 
(A.22) 



In practice, we use the daisy-corrected m 2 2 defined by ( |3.18| ) instead of m\ 2 . 

For the stop with = 0, we should treat the finite-temperature corrections from the 
heavy and light mass eigenstates separately. When mi ^ 0, the stop contributions are 
given by the integrals below. 



Apn 2 3 



A t -A 5 
A t -A 6 



-N c yt/j,A t i 



-N c y 2 tf iA t i f ( } A 1 (k)A 2 (k) 

(-) 



Al(k)Al(k), 



(A.23) 
(A.24) 



(-) 



-A^-gl\A\{k)A 2 {k) 



+ 9 ^A 1 {k)Al{k) +y 2 t ^A\{k)Al{k) 



(A.25) 



A t -A 7 = -N c y*iiA t i [ ( ' 



4 



'g\ 



g!)+y 2 t )Al(k)A 2 (k) 



+ + At(k)^(k) + y 2 t AtA{{k)Ai{k) 



where 



A t {k) 



1 



A 2 {k) 



1 



(A.26) 



(A.27) 



k 2 — m~ ' k 2 
The zero-temperature corrections can also be extracted from these integrals by use of the 
formulas given below. For m~ = 0, we need the infrared cutoff to regularize the zero- 
temperature integrals to obtain the results ( |3.39|) - ( 3.42j) . The finite-temperature Feyn- 
man integrals can be divided into the zero-temperature ones and the finite-temperature 
correction to them, which are usually expressed in terms of one-dimensional integrals. 
The integral appearing in the corrections to m 2 is simplified by the formula 



(T) 



A 1 {k)A 2 {k) 



d 4 k E 1 _L/-(=F)fl^i m2 \ 

(2tt) 4 (k 2 + m\)(k 2 + m 2 2 ) 2^ h \T ' T J 



L(mi,m 2 ) ± / 2 V " 



T ' T 



(A.28) 
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where the zero-temperature part is explicitly evaluated to yield 



L(m 1 ,m 2 ) 



16n 2 



ni 



log- 



m{ — m 2 . M r 2 en m{ — m 2 . M r 2 



> 0, 



(A.29) 



which is renormalized by the DR-scheme just as 
part is given by 



f% (a, b) in the finite-temperature 



ft\a,b) 



for a 7^ b and 



a 2 - b 2 



00 / x' 
dx[ 



x 



Vx 2 + a 2 e ^+^ T 1 v / ?TFe^^l 



1 /•«> 1 



1 



(A.30) 
(A.31) 



2 Jo Vx 2 + a 2 e ^+^ T 1 ' 
/| (a, &) is positive for any (a, 6). 

The integrals in the corrections to Xq^ are calculated by use of the following integral: 

■■ • • " (A.32) 

(A.33) 



where 



H(a) 



with H(l) = -1/2 and 



f!r\a,b) 
for a 7^ 6 and 

(a, b) is negative for any (a, 6). 



a — 1 \a — 1 

dx 



log a — 1 



1 



2{a 2 -b 2 )J ^2T^2 e Vx^ T i a 2 -6 2 



/ 2 W (o,6), (A.34) 



<ix 



3 \/ z 2 +a 2 



+ 



X^ 



+ a 2 ) 3 / 2 ev^? zp 1 x 2 + a 2 ( e ^+^ T I) 2 



(A.35) 



The integral in the corrections to A5 is reduced to 



(T) 



where 



1 ,K 4 U 



1 A^)( m l m 2 
J4 TFTi TFT 



{a - l) 2 Va - 1 
with K{1) = 1/6 and 

1 } 1 /"°° / I 



167r 2 mf V m i/ 27r 2 T 4,/4 VT'T 
a fa + 1 

log a — 2 



(A.36) 
(A.37) 
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2(a 2 -b 2 ) 2 Jo Vx 2 + o 2 e^Tl 



Vx 2 + 62 e^Tl ; (a 2 -6 2 ) 2 



/ 2 W (a,6), (A.38) 
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for a 7^ b and 



fi (a, a) = — dx 
16 Jo 



{x 2 + a 2 ) 5 / 2 e^+^ '=Fl L(a; 2 + a 2 ) 2 



1 



3(x 2 + a 2 ) 3 / 2 J( e v / S^ T l)2 



2 e 2V2: 2 +a 2 

+ 3(x 2 + a 2 )3/ 2 ( e ^+a 2 T i) 3 j- ( A - 39 ) 

f^ T \a,b) is positive for any (a, 6). 

References 

[1] For a review see, A. Cohen, D. Kaplan and A. Nelson, Ann. Rev. Nucl. Part. Sci. 43 

(1993) 27. 

K. Funakubo,Prog. Theor. Phys. 96 (1996) 475. 

[2] M. Carena, M. Quiros, C.E.M. Wagner, Phys. Lett. B380 (1996) 81. 

[3] D. Delepine, J.-M. Gerard, R. Gonzalez Felipe, J. Weyers, Phys. Lett. B386 (1996) 
183. 

[4] P. Huet and A. Nelson, Phys. Rev. D53 (1996) 4578. 

M. Aoki, N .Oshimo and A. Sugamoto, Prog. Theor. Phys. 98 (1997) 1179. 
M. P. Worah, Phys. Rev. Lett. 79 (1997) 3810. 

[5] A. Nelson, D. Kaplan and A. Cohen, Nucl. Phys. B373 (1992) 453. 

[6] K. Funakubo, A. Kakuto, S. Otsuki, K. Takenaga and F. Toyoda, Phys. Rev. D50 

(1994) 1105. 

[7] K. Funakubo, A. Kakuto, S. Otsuki, K. Takenaga and F. Toyoda, Prog. Theor. Phys. 
94 (1995) 845. 

[8] K. Funakubo, A. Kakuto, S. Otsuki, K. Takenaga and F. Toyoda, Prog. Theor. Phys. 
93 (1995) 1067. 

[9] K. Funakubo, A. Kakuto, S. Otsuki and F. Toyoda, Prog. Theor. Phys. 96 (1996) 
771. 

[10] K. Funakubo, A. Kakuto, S. Otsuki, F. Toyoda, Prog. Theor. Phys. 98 (1997) 427. 

[11] D. Comelli, M. Pietroni and A. Riotto, Nucl. Phys. B412 (1994) 441. 

[12] J.R. Espinosa, J.M. Moreno and M. Quiros, Phys. Lett. B324 (1994) 181. 

[13] N. Maekawa, Phys. Lett. B282 (1992) 387; Nucl. Phys. Suppl. 37A (1994) 191. 

[14] L. Dolan and R. Jackiw, Phys. Rev. D9 (1974) 3320. 

[15] S. Coleman and E. Weinberg, Phys. Rev. D7 (1973) 1888. 



26 



